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Abstract 

We analyze the limitations imposed by photon counting statistics 
on extracting useful information about MACHOs from Earth-based 
parallax observations of microlensing events. We find that if one or 
more large (say 2.5 m) telescopes are dedicated to observing a MA- 
CHO event for several nights near maximum amplification, then it 
is possible, in principle, to measure the velocity of the MACHO well 
enough to distinguish between disk and halo populations for events 
with wi m 2 > lday -1 , where oj~ 1 denotes the timescale of the event 
and A m denotes its maximum amplification. Thus, if it turns out 
to be possible to reduce all other sources of error to below that of 
photon statistics, it may be feasible to do useful Earth-based parallax 
observations for high amplification events. 

There are now a number of experimental efforts underway to detect mi- 
crolensing by massive compact halo objects (MACHOs). |1|, 0, [| To date 
there have been over 50 events reported. Events are now being discovered 
in realtime, ||, || and telescopes can be trained on events in the process of 
enfolding. 

It is of considerable interest to determine the masses and velocities of 
the MACHOs, since from the masses one could determine their contribution 
to the dark matter in our galaxy, and from their velocities one could de- 
termine whether they belong to a disk or halo population. Unfortunately 
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these parameters are not, in principle, determinable by observations from a 
single telescope undergoing inertial motion. This can be seen as follows: the 
amplification of a lensed star (which is all that can be measured) is given by 

a( \ 1)2 + 2 m 
MU) = + 4)V2 ' (!) 

with u = d/r e , where d is the perpendicular distance of the MACHO from the 
line joining the observer to the star. Here r e is the Einstein radius, defined 
by 

n = ^ Gm f L ~*\ ( 2 ) 

where m is the mass of the lensing object, L the distance from the observer 
to the star, and x the distance from the observer to the lens. For inertial 
motions of the star, MACHO, and observer, the time variation of u is given 
by 

u(t) = vW^(t-to) 2 , (3) 

with to the central time, u> = v/r e , and (3 = b/r e , where v is the velocity of 
the MACHO with respect to the line between the observer and the star, and 
b is its impact parameter (so that uj _L (3). Consequently, the light curve of 
a MACHO event is completely characterized by the three parameters (A m , 
uj, t ), where A m is the maximum amplification, given in terms of f3 by, 

P 2 + 2 

Am = f3(^ + 4)1/2' ( 4 ) 

It follows that the only information about v and m that can be extracted 
directly from the light curve is the value of the single parameter uj = v/r e . 

The situation improves considerably, however, if an additional telescope 
(widely separated from the first) also observes the event. 0, H (Alterna- 
tively, similar improvement can occur if a single telescope undergoes signifi- 
cant deviations from straight line motion over the timescale of the event. [0, 
||) As a result of parallax effects, the second telescope will see a similar light 
curve, with the same value of uj, but with slightly different values of A m 
and to- As we shall see explicitly below, from this information one can, in 
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principle, [] determine the "reduced Einstein radius" 



r ~ e = Jl^x)^ (5) 

The magnitude of the "reduced transverse velocity" (v = r^z^j v) is then 
given by 

v = — -r e u = r e u. (6) 

[L — x) 

The parameters f e and v yield considerable further information about m and 
v, although one must still know or estimate the value of x/L from other 
arguments to determine m and v exactly. 

For separations as large as 1 Au there should be no difficulty in measuring 
the above parallax effect. Indeed, parallax effects due to the Earth's orbit 
about the Sun have recently been observed in an event lasting over two 
hundred days (the longest event yet detected). || However, in order to make 
use of the Earth's orbital motion, the lensing event must span a fair portion 
of the earth's orbit. Most events observed thus far have considerably shorter 
timescales (< 100 days), and therefore do not lend themselves to this type of 
measurement. With an additional telescope placed in solar orbit, however, 
it should be possible to measure parallax effects for short timescale events 
without difficulty, as has been proposed by Gould. [RL 

Unless one undertakes the major project of putting a telescope into orbit, 
the effective baseline for a parallax observation of a short timescale event 
(achieved either via an additional telescope or via the rotation of the earth) 
is limited to about an earth radius. It has generally been assumed that the 
very small changes to the light curve associated with such a short baseline 
would be insufficient to allow any useful parallax information to be extracted. 
However, the only truly fundamental limitation on the ability to determine 
the parameters of a light curve arises from "photon statistics", i.e., the ran- 
dom y/N fluctuations in the number of photons reaching a telescope. In this 
paper, we investigate the limitations imposed by photon statistics on parallax 

1 As analyzed in there is a 2-fold degeneracy in the determination of f e , as it 
is possible for the trajectory of the line between the lensed star and MACHO to pass 
between the telescopes. However, the probability of this happening for the short baseline 
observations of interest here is completely negligible, and thus will be ignored in our 
analysis. (In any case, the degeneracy can be resolved by using three telescopes or two 
telescopes undergoing non-inertial motion. pQ|) 
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observations using Earth-based telescopes. Our conclusion is that although 
these limitations are formidable, it is possible in principle to extract useful 
information for high amplification events. Specifically, with one or more large 
(say 2.5 m) telescopes observing for several full nights, it should be possible 
to determine v well enough to distinguish between disk and halo populations 
for events with uM m 2 greater than about lday -1 . Thus, if all other sources 
of error can be adequately controlled, it may be possible to obtain some new, 
valuable information about MACHOs using Earth-based astrometry. 

The limitations imposed by photon statistics can be estimated by analyz- 
ing the following two idealized problems for a single inertial telescope which 
makes perfect measurements of the photon flux reaching it: (i) Suppose the 
correct values of the parameters uo and A m of the light curve for a MACHO 
event are given. What error would photon statistics be expected to produce 
in the determination of to from observations of the light curve? (ii) Suppose, 
similarly, that u and to are given. What error would photon statistics be 
expected to produce in the determination of A m ? 

To answer the first question, we imagine "binning" the photons which 
arrive at the telescope into intervals of time T. The expected photon flux in 
the bin at time t is then N = foA(t)T, where fo denotes the average photon 
flux from the unamplified star and A(t) denotes the amplification produced 
by the MACHO at time t. In principle the central time, to, could be obtained 
from a measurement of N in a single bin (say the ith bin) by inverting this 
equation, using eqs. |l] and || together with the given values of u and A m . 



However, photon statistics produces an error (AiV)j = y/N{ = \J JqA^T. 

This causes an error, (Ato)«, in the determination of to from a single-bin 
measurement given by 



Since the photon statistics errors for the different bins are uncorrelated, the 
total error resulting from making measurements of this type for each bin on 
the entire light curve is given by 





-i -l 
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(Ai °» 2 = I? (I) 
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dt 



(8) 
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where in the final expression we have taken the limit T — > 0. By similar 
reasoning, the answer to the second question is that the error due to photon 
statistics in measuring the maximum amplification, for given u and t ? is 

-l 



(AA T 



(9) 



The integrands of eqs. [8] and [5] are plotted in Fig. [I], for the light curve 
corresponding to the A = 7 event observed by the MACHO group. [|IJ 

Equations [] and |9] simplify considerably for events where A m is large - as 
will be the only case of interest here - since then eq. [I] is well approximated 
by A(u) = 1/u over the relevant portion of the light curve. Substituting 
A(u) = 1/u into eqs. [8] and |9] and performing the integrals, we obtain, for 
large amplification events, 

(Ai °» 2 = v^aj (10) 

and 

(AA m ) 2 = *^=L. (11) 

In the large amplification approximation we find that the the integrand of 
eq. H has a double peak at times tp given by 



'2/3 

t P = t ± ^— , (12) 

LOA m 

where the width of each peak is of order l/uA m . The integrand of eq. |9| has 
a single peak at t — to, with a width also of order l/ujA m . 

Equations |8] and [], and their approximations |10| and [11], provide lower 



bounds on the errors arising in measurements of to an d A m by a single tele- 
scope. Clearly, useful information about the MACHO can be extracted from 
parallax observations by two telescopes only if the errors in measuring to 
and/or A m by the individual telescopes are small compared with the differ- 
ences, Sto and 5A m , occurring in the values of these parameters between the 
two telescopes. If we have two (inertial) telescopes separated by a transverse 
distance D, then the shift in the central time between the telescopes is given 
by (see footnote 1) 

D cos 9 

St = — 13 

ijjr P 
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where 9 is the magnitude of the angle between v and the line joining the 
telescopes. Similarly, the shift in ft is given by 



5ft 



Dsin6 



(14) 



In the large amplification approximation (where A Tl 

= _D4 a 9A m > 



1 /ft) this yields 



(15) 



Thus, the appropriate "figures of merit" for parallax measurements - as- 
sumed to be limited only by photon statistics - for large amplification events 
are 

'Af n \ 2 3* 



and 



St 
/AA V 



2f ujA m 2 D 2 cos 2 9 
3d 2 



(17) 



\5A m J Af tuAjD 2 S m 2 9' 

Remarkably, apart from the replacement of cos 2 9 by sin 2 9, the right sides 
of eqs. [IB] and [Tj] differ only by a factor of 2. This means that - apart from 
different angular sensitivities - the time delay and amplification differences 
are essentially equally useful for measuring v. 
Solving eqs. IT3 and O for v = uir e , we obtain 



v 2 



DWAJ 



{SA m ) 2 + u 2 A m \6t f 



-i 



(18) 



For the purpose of making an order of magnitude estimate, we shall take the 
errors caused by photon statistics in the determination of to and A m as being 
independent, and given (for each telescope) by eqs. |TU] and [TT| respectively. 
The resulting error in the determination of the magnitude of v from a parallax 
measurement by two inertial telescopes is given by 



Av 

v 



V 

15 



2f uA r 



1/2 



VT+ cos 2 ^. 



(19) 



Quite valuable information would be obtained from a parallax measure- 
ment if it were able to determine v accurately enough to distinguish between 
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a disk and halo population. As a rough criterion for this, we require that 
one be able to make a ~ 50% measurement of v when v ~ 80 km/sec. Let 
us estimate the conditions which must be satisfied by a lensing event in or- 
der for this to occur. We consider a typical lensed star with an unamplified 
magnitude of V=19. For a 2.5 m telescope and an observing bandwidth of 
1000 A, this corresponds to fo = 10 3 photons/sec. For the purpose of making 
an estimate, we assume that two 2.5 m telescopes are dedicated "full time" to 
the task of observing the MACHO microlensing event for, say, 3 or 4 nights 
near maximum amplification. The effective photon flux is reduced in the 
above formulas by a factor of 3 to account for an 8 hour per night observing 
period, and by another factor of about 2 to account for our observing only 
near maximum amplification. (Since the width of the peaks of the integrands 
in eqs. |] and |9] vary as l/uA m , a somewhat larger reduction would have to 
be made for relatively low amplification or long timescale events.) Choosing 
D to be an Earth radius (6, 000 km) and taking cos 2 9 = 1/2, we find that to 
distinguish between disk and halo populations in our two-telescope idealized 
parallax experiment, we need Q 

^C^lday- 1 . (20) 

Thus, for a typical value of u = .05day _1 , we see that the errors resulting 
from photon statistics alone will limit useful measurements to the cases where 
A m > 5. For Earth based measurements, this limit could be improved only 
by use of more (or larger) telescopes. Note that since, for any given MACHO 
mass and velocity, the probability of an event having impact parameter less 



than b is proportional to 6, [1 L] and since for large A m we have A m ~ 1//3, it 
follows that about 20% of all events will have A m > 5 (where an "event" is 
defined by the criterion that (5 < 1). Most importantly, the imposition of a 
selection criterion on A m alone randomly samples the MACHO events, i.e., 
it does not bias the class of all MACHO events with respect to any MACHO 
properties. However, a restriction to events within a range of timescales can 



2 Note that the lower limit for toAm in eq. gC] scales with D as D~ 2 . Hence, if the 
second telescope were placed in solar orbit at a distance of 1 Au from Earth, one would 
decrease this limit by a factor of 5 x 10 8 (although it should be remembered that eq. ^ 
applies only in the high amplification approximation). Thus, even taking into account the 
much smaller area of an orbiting telescope, we see that photon counting statistics would 
not impose any significant limitations on parallax measurements using a telescope in solar 
orbit. 



7 



introduce biases. 

Our analysis above has ignored the rotation of the Earth. In fact, the 
Earth's rotation is helpful in that it enables one, in principle, to do parallax 
measurements with a single telescope (since it provides an effective baseline 
of order an earth radius) - although, undoubtedly, one would want to do the 
measurements with at least two telescopes to improve statistics and reduce 
systematic errors.0 Nevertheless, one might worry that the estimates used in 
deriving eq. [19] could be altered by effects of the Earth's rotation. In addition, 
the estimates of Ato and AA m in eqs. |8] and |9| may have been optimistic in 
that they were derived under the assumption that all the parameters except 



the one in question were known exactly. To check the validity of eq. [19 
we have written a Monte Carlo code in which one or more telescopes were 
placed at various geographic locations on a rotating Earth, each observing 
for 8 hours per night for several nights. The theoretical light curve for each 
telescope was then binned into 0.2 hour "observations", and each observation 
was then randomly altered by the appropriate Gaussian \/N photon statistics 
errors. We then fit the resulting observations with the six free parameters f3, 
u>, to? T e . Repeating this for a large number of trials, we arrived at numerical 
estimates for At , Av4 m , and Av. These numerical estimates were compared 
with the corresponding analytic estimates, obtained by utilizing appropriate 
"effective baselines" and modifying the range of the integrals in eqs. [8] and |9] to 
cover only the period of time during which the telescopes actually observed. 
The agreement between the numerical and analytic estimates was very good 
(within a factor of 2 in essentially all cases tried). Thus, we are confident that 
the formulas given above provide reliable order of magnitude estimates of the 
limitations imposed by photon statistics errors on parallax measurements of 
MACHOs. 

The key issue regarding the feasibility of Earth-based parallax measure- 
ments of high amplification events (A m > 5) is whether all other possible 



3 Note that with two telescopes separated by latitude on a rotating Earth, one has 
effective baselines in both latitude and longitude. Consequently one could, in principle, 
measure both components of v by making use only of the information contained in SA m . 
Since, when suitably normalized, the peak of the (AA m ) 2 integrand tends to be higher and 
sharper than the double peak of the (Aio) 2 integrand (see Fig. 1), this suggests that the 
best strategy for parallax observations might be to have many large telescopes (separated 
in latitude) observe the event for one or two nights near maximum amplification, rather 
than have one or two telescopes observe the event for many nights. 
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sources of error can be reduced to below that of the photon statistics error 
computed above. The challenges to doing this are quite formidable. How- 
ever, it should be noted that some of the conditions present in a microlensing 
event are extremely favorable for doing accurate photometry. In particular, 
the brightness change of the lensed star is completely achromatic, and the 
required brightness measurements are differential both with respect to time 
and with respect to the background stars in the field. Since the lensing 
events are observed in crowded star fields, there should be no difficulty find- 
ing an appropriate "mix" of nearby stars to match the color of the lensed 
star. Thus, it should be possible to eliminate the effects of time dependent 
and color dependent changes in the absorption by the atmosphere. Unfor- 
tunately, crowded fields also exacerbate what is probably the most serious 
source of error: the time dependence of seeing conditions, which causes a 
spread in the stellar images. With the high magnifications attainable by 
large telescopes (which would have to be used for the parallax measurements 
in any case), it should be possible to accurately measure, and correct for, 
the effects of seeing conditions on stars in a field. However, we must leave a 
detailed analysis of sources of error and other practical limitations to those 
more experienced than we in such matters. 
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9220644 to the University of Chicago. D.H. acknowledges support from an 
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Figure 1: The integrands of eqs. |8] and[^, plotted for the first event reported 
by the MACHO group (A m = 6.86, 1/u = 16.95). The (At ) 2 integrand 
has been multiplied by u; _2 y4 m _4 in order to make it dimensionless; this also 
weights it appropriately relative to the (AA m ) 2 integrand with respect to 
utilizing to and A m to obtain information about the velocity of the MACHO 
(see eq. [18] below). From the figure, we see that essentially all of the infor- 
mation concerning A m can be extracted from observations made during a 
period of a few days near maximum amplification, whereas the information 
concerning to is peaked (also with a spread of a few days) two days before 
and after maximum amplification. 
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